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Abstract 

A connection between the Yang-Baxter relation for maps and the multi-dimensional consis- 
tency property of integrable equations on quad-graphs is investigated. The approach is based 
on the symmetry analysis of the corresponding equations. It is shown that the Yang-Baxter 
variables can be chosen as invariants of the multi-parameter symmetry groups of the equations. 
We use the classification results by Adler, Bobenko and Suris to demonstrate this method. Some 
new examples of Yang-Baxter maps are derived in this way from multi- field integrable equations. 



1 Introduction 



The quantum Yang-Baxter (YB) equation has its origins in the theory of solvable models in statis- 
tical mechanics 121 and the quantum inverse scattering method ^ . The fact that this equation 
has also found many applications in representation theory, the construction of invariants in knot 
theory and that it lies at the foundation of quantum groups, gives to the quantum YB equation 
a prominent position among the basic equations in mathematical physics, see e.g. |5, (Hj and 
references therein. 

In its original form, the quantum YB equation is a relation for a linear operator R : V ^ 
V where F is a vector space. The relation has the form 



in End{V (>i)V (SiV), where R^^ is meant as the identity in the second factor of the tensor product 
y (g) y (g) y and as R in the first and third factors, and analogously for R^^, R?^ . Supposing that 
X is any set, the maps R from the Cartesian product X X into itself, which satisfy the relation 

of composite maps are called set theoretic solutions of the quantum YB relation. The study of 
set theoretic solutions of the quantum YB equation was originally suggested by Drinfeld [S] (see 
also earlier work by Sklyanin jjj, where the first interesting example of such solution was found) 
and since then they have attracted the interest of many researchers. 

More recently, a general theory on the set theoretic solutions to the YB relation was developed 
in [Hj and the notion of transfer maps, which can be considered as the dynamical analogues of the 
monodromy and transfer matrices in the theory of solvable models in statistical mechanics, was 
introduced in . In many interesting examples of YB maps ^ , such as maps arising from geometric 
crystals jlU| , the set X has the structure of an algebraic variety and is a birational isomorphism. 
The case of CP^ X CP^ has been recently discussed in in relation with the classification of the 
so-called quadrirational maps. 

In this paper we investigate the relation between the YB property for maps and the multi- 
dimensional consistency condition for equations on quad-graphs, which is now commonly accepted 
as a definition of integrability for such equations (see jl2| I14j). Although the link between 
these two notions was known before (see e.g. concluding remarks in it was never explored 
systematically. 

Our approach is based on the symmetry analysis of integrable equations on quad-graphs. The 
main idea is that the YB variables are suitable invariants of their symmetry groups. A good 
example is the discrete potential Korteweg - de Vries equation (dpKdV) ^1 El 



(see notation on Fig. It is clearly invariant under the translation / — > /-|- const. The invariants 

^Adopting the terminology in j^, set theoretic solutions to the YB equation will be referred in the following simply 
as YB maps. 
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X = fl-f, y = /l,2-/l, U = fl,2-f2, V = f2-f, (3) 

satisfy the relation 

X + y = u + V. (4) 
and equation ^ is written in terms of them as 

{x + y){x - v) = ai - a2. (5) 

This aUows to express u,v as functions of x,y, which leads to the following YB map 

u = y^ ■ , v = x ■ , (6) 

x+y x+y 

known as the Adler map |^. Note that the YB variables x, y, u, v, are attached to the edges of the 

lattice. The fact that the corresponding map satisfies the YB property follows directly from the 

3D consistency property of dpKdV (see Fig. This construction works for integrable equations 

on quad- graphs with one-parameter symmetry group. 

One of the main findings of our paper is that this idea works for multi-parameter symmetry 

groups if one considers the extension of the equation on a multi-dimensional lattice. In that case 

the edges are replaced by higher dimensional faces. We show that in such a way one can derive 

from the same discrete potential KdV the following YB map 

u = yQ. „ = Q, ('-7-2) + fe-7.)^ + 72(7l-l):^./ ^ ( 

(1 -7i) + (71 -72) 2/ + 71(72 -^)xy 

We will call it Harrison map since it is closely related to the superposition formula of the Backlund 
transformation for the Ernst equation in general relativity introduced by Harrison ^Hl- After the 
change of variables xi-^l/x, v 1/v, yi— >y/72, ui-^n/7iit coincides with the Fj quadrirational 
map in which corresponds to the most general case of two conies. Note that the most degenerate 
case Fy in the classification of is simply related to the Adler map. 

The plan of the paper is the following. We start in section[21with the discussion of 3D consistency 
property for the equations on quad-graphs. As the examples we choose three equations from the 
classification list in jl4j . By considering the invariants of their one-parameter symmetry groups 
we derive all five types of the quadrirational maps from Next, in section |21 we show how 

this symmetry method can be generalized in the case where the lattice equation admits a multi- 
parameter symmetry group. This is demonstrated on the example of the lattice KdV equation by 
extending it to a three dimensional cube and using the invariants of a two-parameter symmetry 
group as YB variables. Finally, we show how the Harrison map can be retrieved from the lattice 
KdV equation by exploiting its full three-parameter symmetry group and the consistency property 
on a four dimensional cube. 

In section |3] we show that the same idea works equally well for multi-field integrable lattice 
equations and we derive certain examples of multi-component YB maps from lattice equations 
in the Boussinesq family, vector Calapso equation and its specialization to an integrable discrete 
version of the 0{n + 2) nonlinear cr-model, introduced recently by Schief [20]. The paper concludes 
with perspectives where we address some questions for future study. 
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2 Multi-dimensional consistency and YB maps 



2.1 Equations on quad-graphs and the 3D consistency property 

Central to our considerations are integrable discrete equations on quad- graphs, which are specific 
equations associated to planar graphs with elementary quadrilaterals faces. In the simplest case 
one has complex fields / : — > C assigned on the vertices at sites (ni, 712) and two complex lattice 
parameters ai, 02 assigned on the edges of an elementary square being equal on opposite edges (see 
Fig. ^ . The basic building block of such equations consists of a relation of the form 

^(/,/i,/2,/i,2;ai,a2) = 0, (8) 

between the values of four fields residing on the vertices of each elementary quadrilateral for which 
we use the shorthand notation: 

/:=/(ni,n2), /i := /(ni + 1, 7^2), /2 := /(ni, n2 + 1), /i,2 := /(ni + 1, n2 + 1) . (9) 



fl,2 



a2 



a2 



fi 



Figure 1: An elementary quadrilateral. 

Integrable discrete equations of the above type © are listed in a recent classification Jl], where 
the 3D consistency property (see below) and some additional conditions were imposed. Prom that 
list we consider the following equations 



£1 
£2 
£3 



(/i,2-/)(/i-/2)-ai + a2 = 0, 

Mf fi + f2 fi,2) - 02(7 /2 + fi fi,2) + 5(ai' - 02') = , 

(1 - a2^)(/i - ai/)(/2 - ai/1,2) - (1 - ai^)(/2 - a2/)(/i - 02/1,2) 



0. 



(10) 

(11) 

(12) 



Equation £1 is already mentioned as the dpKdV equation. Equation £2 with (5 = is the 
modified discrete KdV or Hirota equation 15 . If 7^ we may always assume that 5 = 1 using an 
appropriate gauge. Equation £^ corresponds to the equation labeled as Q?>s=q in the classification 
of IE]- It is contained in the 4-parameter family of the equations derived earlier in |21j . which 
contains also discrete versions of potential KdV, modified KdV and Schwarzian KdV (see |22] for 
a more recent discussion). 

The integrability of such equations can be defined using the three dimensional consistency 
property. This means that the overdetermined system consisting of the difference equations 



£if,fiJ'j,fi,j;oii,aj) = 0. 



1 < i < j < 3 , 



(13) 
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and their shifted versions, is consistent on the three-dimensional lattice . In practice, this 
property is verified as follows |12| I13[ I14j . Consider an elementary initial value problem on the 
three-dimensional cube with initial data assigned on four vertices, not all of them lying on the 
same face. One such initial configuration is depicted in Fig. |2{a) with initial values /, /j, 1 < i < 3. 
Using equations ()13() on the three faces adjacent to the vertex with value /, we determine uniquely 
the values /jj-, 1 < i < j < 3, in terms of the initial data. Then using shifted versions of (|13|) on 
each of the remaining three faces, we evaluate /i,2,3 in three different ways. Consistency means 
that one obtains the same value for /i,2,3 in terms of the initial data /, /i, 1 < i < 3 (independent 
of the way we choose to evaluate it). For the dpKdV equation (jJJ this value is 




Figure 2: Elementary initial value problems on the cube 



^ ^ _ («! - «2)/l /2 + (Q3 - /3 + («2 - Q3)/2 ^?, ^^^^ 

(a2 - ai)/3 + (ai - a3)/2 + (03 - a2)/i 

Note that the right hand side of equation H14|) is invariant under any permutation of the indices 
(1,2,3) which label the field variables and the associated lattice parameters. 

Another initial data configuration, which is best adapted to the YB property that we consider, 
is depicted in Fig. Efb). A third possible initial configuration is to give the values /, /i, /2, /i,2,3- 
The latter two configurations are equivalent to the first one by using the equation on one of the 
faces. For example, by using the front face equation we can exchange the value /i_3 to /s in the set 
of initial data. 

Using the fact that dpKdV equation possesses the 3D consistency property, one can show in a 
similar manner that it can be consistently imposed on each 2-dimensional face of a 4-cube. Since 
we are going to use this property later on we describe explicitly its derivation. 

For given initial values /, i = 1,2,3,4, we determine the shifted values of the fields involving 
any two different directions, using the equations 

iki - f){fi- fj) = oii- aj , (15) 

1 ^ ^ < J ^ 4 (see Fig. Successively, since dpKdV is 3-dimensional consistent, we determine 
the values fijk, l<i<j<k<4. Then the value /i,2,3,4 can be found in six different ways, 
using the dpKdV equations on the six 2~dimensional facets containing the vertex where the value 
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/i,2,3,4 is assigned. This vertex is contained also in four cubes, each one of them containing three of 
the six facets, and the incidence relations are such that taking into account the three dimensional 
consistency on each of the four cubes one proves that the value /i,2,3,4 is uniquely determined in 
terms of initial data. By direct calculations also we find that this value is independent of the way 
that we used to calculate it, and equals 



fl,i,j,k 



a {aiaifij + ajakfj,k){fi - fi){fj - fk) 



cr {aiUi + ajak){fi - fi){fj - fk) 

ijk 



(16) 



where cr denotes the cyclic sum over the subscripts {i,j,k) = (2,3,4), (4,2,3), (3,4,2). It can be 

ijk 

easily checked that /i,2,3,4 given by (|16|) . remains invariant under any permutation of the indices 
(1,2,3,4), thus dpKdV is four dimensional consistent. 




1,3,4 



/2,3,4 /l, 2,3,4 

Figure 3: Discrete potential KdV in 



2.2 YB relation and 3D consistency property 

Let X be any set and R a map of X x X into itself. Let R'^ : X" ^ X", where X" = X x X x . . . x X, 
denotes the map which acts as R on the i and j factors and as the identity on the others. More 
explicitly, let us write R{x,y), x,y S X, as 

R{x,y) = {f{x,y),g{x,y)). (17) 

Then, for n > 2 and 1 < i, j < n, i ^ j we define 

R'^ {x\ X'^ , . . . x"") = { (18) 

{x^, . . . , x^~'^,g{x'^,x^),x^^^, . . . , x*"-^, /(x*, x^),x'^~^^, . . . x") i > j ■ 
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In particular, for n = 2 we find that R^^ = R and R'^^{x,y) = (^g{y,x), f{y,x)). The latter map 
can be written as a composition of maps as follows 



R^^ = PRP. 



(19) 



where P is the permutation map, i.e. P{x,y) = {y,x). 

A map R is called a YB map if it satisfies the YB relation , regarded as an equality of maps 
of X X X X X into itself. If in addition the relation R^^ i? = Id holds, then R is called reversible YB 
map. 

In a more general setting we may consider a whole family of YB maps parametrized by continu- 
ous parameters rather than a single map. The YB relation then takes the parameter-dependent 
form 

i?23(a2, as) i?^^(ai, as) -R^^(ai, as) = i?^^(ai, aa) i?^^(ai, ag) R^^{a2, 03) , (20) 
and the reversibility condition becomes 

i?2^(a2,ai)i?(ai,a2) = Id. (21) 

The relation between YB maps and integrable equations on quad-graphs can be demonstrated 
in the example of the discrete potential KdV equation £1. As we have already shown in the 
Introduction, by considering the differences of the values of the fields assigned on two adjacent 
vertices ©, we arrive at the Adler map (jH)). There is a different combination for the variables 



{v,a2) 




(y,a2) 



(x,ai) 



Figure 4: An oriented quadrilateral for the map R{x,y) = {u,v). 
assigned on the edges of the square, namely 

X = ffl, y = /l/l,2, U = f2fl,2, V = ff2. 

From the above relations (1221) we deduce that 



xu = yv . 

Moreover, dpKdV can be also written in terms of the variables ()22() as follows 

y + v — X — u = ai — 02- 



(22) 



(23) 



(24) 
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Figure 5: Three dimensional representation of the YB relation 

Solving equations (|23|) and (|24j) for (u, v) we get the following map 

M = y 1 H , V = X . (25) 

V x-y J V x-y J 

The maps automatically satisfy the parameter dependent YB relation (|^n)l . Indeed, 

it is easily shown that the consistency property for a configuration of initial data on the vertices 

of a cube, as depicted in Fig. \^ is equivalent to that with initial values f , fi, fi.2, f 1,2,3- These 

initial data correspond to the values {x,y,z) on the edges (Fig. ISJ. The 3D consistency property 

guarantees that the composite maps 

(a) : (x, y, z) — > {x2, z) — > (x2,3, V-i, z-2) — > {x2,3,y-i,3, z-1,-2) (26) 

(b) : (x, y, z) — > {x, ys, z_2) — > {x^, ys, -2) — > (^2,3, y-1,3, Z-1-2) (27) 

appearing in equation (^OJ, applied on {x,y,z) give identical values for (x2,3, y-i, 3, -2-1,-2)- 

Analysing these two examples one notices that the variables x, y, u, v, which we call YB variables, 
are invariants of certain symmetry groups of the relevant lattice equation. Now we are going to 
show that this symmetry method can be applied in more general situations as well. 



2.3 Lattice invariants of symmetry groups and YB variables 

Let us first recall the basic notions of Lie symmetry methods applied to lattice equations of the 
form (jHl). With minor modifications these are in accordance with the symmetry methods applied 
to algebraic or differential equations (see e.g. ^1 for an extensive study on the subject). 

Consider a lattice equation of the form Q involving one field / : ^ C (CP^). Let G be a 
one-parameter group of transformations acting on the domain of the dependent variables, 

G:/^$(ni,n2,/;e), egC. (28) 

The prolongation of the group action on the lattice jet space J with coordinates (/, /i, /2, /i,2) is 
specified by 

G : (/,/i,/2,/i,2) ^ (^>(ni,n2,/;e),$i(ni+l,n2,/i;e),^>(ni,n2+l,/2;e),'^'(ni+l,n2+l,/i,2;e)) • 

(29) 
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The infinitesimal generator of the group action of G on / is the vector field 

(30) 



Qini,n2,f)df , where Q{ni,n2, f) = ^'^{ni,n2, f; e] 



e=0 

There is a one-to-one correspondence between connected groups of transformations and their as- 
sociated infinitesimal generators since the group action is reconstructed by the flow of the vector 
field V by exponentiation 

^(f^i,"2,/;e) = exp(ev)/ . (31) 

The prolongation of the infinitesimal action of G given by (|29j) , is generated by the prolonged vector 
field 

v = Qdf + Qidf,+ Q2 df, + Qi,2 df,, , (32) 

where subscripts denote Qi = Q{ni + 1,^21/1)1 Qi,2 = Qin-i + l,n2 + l,/i,2)i and so on. 

The transformation G is a symmetry of the lattice equation ©, if it transforms any solution 
of Q to another solution of the same equation. Equivalently, G is a symmetry of equation Q, if 
the equation is not affected by the transformation ()29|). The infinitesimal criterion for G to be a 
symmetry of equation is 

v(£'(/,/i,/2,/i,2;ai,a2)) = 0, (33) 

whenever equation (jSJ holds. 

A function / : J — > C is a lattice invariant of the transformation group G, if / is not affected 
under the action of G. The infinitesimal invariance condition for the lattice invariants is 

v(/) = 0. (34) 

Once we have determined a symmetry generator v of the lattice equation Q, the corresponding 
lattice invariants can be found from the solution of the first order partial differential equation (|34|) , 
by using the method of characteristics. From the corresponding system of ordinary differential 
equations we may easily obtain the general solution, since it consists of equations with separated 
variables. We assign now to the edges of an elementary quadrilateral the following YB variables 
(Fig.® 

X = /(/,/i), y = /(/i,/i,2), U = /(/2,/l,2), v = Iif,f2), (35) 

where / is an invariant depending on two neighboring values of /. Since G is a symmetry of the 
lattice equation, the latter can be written in terms of these variables: 

X'(x, y, u, v; ai , 02) = . (36) 

This can be done in different ways since the variables (|35|1 are not independent; there exists a 
relation among them 

J^{x,y,u,v;ai,a2) = , (37) 

following from the fact that the space of G-orbits is three-dimensional. 

Solving the system of equations (|S7)) for u, v in terms of x, y and assuming that the solution 
is unique, we obtain a map R{x, y) = (n, v). 
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Proposition 2.1. // the discrete equation 8 satisfies the 3D consistency property, then the map 
R{x,y) = {u,v), which relates the lattice invariants icl5\) . satisfies the YB relation. 

The proof follows from Fig. [21 Similar considerations hold for multi- field lattice equations, 
which give rise to multi-component YB maps (see section^. 

Example 2.2. Consider the dpKdV equation Two infinitesimal symmetry generators of the 
latter equation are 

VI =9/, v2 = (-ir+"^/9/, (38) 
(see [21] )■ They generate the symmetry transformations 

G':f^f + ei, :/^/exp(e2(-ir+"^) , £1,62 gC, (39) 

respectively. The lattice invariants assigned on the edges of a square for each one of the above 
symmetry transformations are the variables © and (|22j) . respectively. 

The consideration of the remaining two equations from the list ()1U() - H12|) . with the corresponding 
symmetry generators 

^£, = {-ir+-^fdj, ^£, = fdf, (40) 

leads to the results summarized in the following Tabled 

Table 1: Yang-Baxter maps arising from equations £2 and £3. 
Yang-Baxter Functional relation 

Yang-Baxter map 

variables and lattice equation 



3 



X 



f fi/ai y=y- 7i(a^ + <^) -72(y + ^) 

£2 _r r I 71 (x + u + ()) = 72 (y + w + d) 

^ - -^^'^/"^ where 7. = „ - ^ 7i (x + ^) - 72 (2/ + <5) 

£ £ I ILL y 

v = J h/a2 72 x-y 



^ — fl/i'^lf) XV = UV r\ 

y u = yQ, v = xQ 

y = /i,2/(a2/i) 1 - x"^ _ 1 - 71 1 - 72 y 

u = /i,2/(ai/2) l-v-^ 1 - 72 1 - 71 « 



V = f 2/(012!) 



where 7i = ^ (i-7i)+(7i-72) j/+7i(72-i) 



(i-72)+(72 -7i) ^+72(71-1) 



These maps are simply related to quadrilateral maps -F/ — i*/// from the Adler-Bobenko-Suris 
list of jllj . Namely, setting 5 = in the YB map constructed from £'2, we retrieve the map labeled 
as Fill map in The case 5 = —1 corresponds to the Fu map. Finally, the YB map constructed 



10 



from 82, under the transformation x ^ 1/x, v ^ V ^ V I li-, u/^i turns into 

PC n n (1 - 72)x + 72 - 71 + (71 - 1)?/ /..s 
u = 7iyQ, v = -i2xQ, Q = — — — ■ — , (41 

72(1 - + (71 - 72)a;y + 71(72 - ijy 
which corresponds to the Fi map in the classification in ill . As we have aheady mentioned it is 
closely related to the superposition formula of the Backlund transformation for the Ernst equation 
in general relativity introduced by Harrison jl9j . c.f. [25| . 

The remaining maps Fiy and Fy from are related in a simple way to the maps © and 
(^5]) derived from dpKdV already in section 12.21 

Thus all 5 types of quadrilateral maps from the Adler-Bobenko-Suris classification are 
equivalent to the YB maps coming from the integrable equations on quad-graphs. 

Remark 2.3. We should mention that the equivalence of the quadrilateral maps considered in jllj 
allows independent change of variables x, y, u, v and therefore does not respect the YB property, 
which is preserved in general only under the diagonal action of the Mobius group. In particular, 
quadrilateral maps in general do not satisfy the YB relation (contrary to what one might conclude 
from ) . However sometimes two YB maps are related by non-diagonal action (see the example of 
Harrison map and Fy above). The question how many such pairs exist needs further investigation 
(see the discussion of this in J6j ) . 



3 Multi-parameter symmetry groups and multi-dimensional con- 
sistency 

The purpose of this section is to show that the symmetry method described in the preceding section 
works equally well for multi-parameter symmetry groups. The idea is to consider the extension of 
the equation into many dimensions using the SD-consistency property and then prescribe the YB 
variables not to the edges but, for example, to higher dimensional faces. 

We demonstrate how the method works in the example of dpKdV equation, which is invariant 
under the three-parameter symmetry group G with infinitesimal generators 

VI = 5/, v2 = (-ir+"^/5y, V3 = (-ir+"^a/. (42) 

Their commutators are 

[V1,V2]=V3, [v2,V3] = -Vi, [vi,V3]=0, (43) 

from which it is immediately seen that {vi,V2,V3} span a solvable Lie algebra. This algebra is 
actually isomorphic to the Lie algebra of the group /so(M^'^) of isometrics of Minkowski plane, so 
G can be considered as the connected component of identity of this group. 

3.1 Consistency of dpKdV around a 3-cube and the Fm map 

Consider the dpKdV equation imposed on each face of an elementary cube (see Fig. \^ and the 
abelian subgroup H of the full symmetry group G generated by {vi,V3} (translations of the 
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Minkowski plane). These two symmetries can be extended to the corresponding system on 1? ^ 
e.g. V3 = (-l)"i+"2+"39^ . 

We consider now the following invariants of the subgroup H 

X = fl- fs, y = fl,2- fl,3 , U = /i,2 - /2,3 , V = f2- fs, (44) 

assigned on four faces of the 3-cube. Using the dpKdV equations 

ai-a2 „ , ai - 03 , 012 - 
Ji,2 - J — —, J- , Ji,3 - J — —f T , /2,3 - / — 7- J V^bj 

J I — J2 Jl — J3 J2 — J3 

we easily find that the invariants (|44j) are related by 

^1 ^2 .... 
uv = xy, u = y , (46) 

X V 

where Pi = ai — a^, P2 = 02 — a^. Solving the system for {u,v) in terms of {x,y) we obtain 
the map 

u = yP, v=xP-\ P=^^^, (47) 

which satisfies the YB relation as it can be checked by direct calculations. This fact is also related 
to the higher dimensional consistency of dpKdV on as it is explained in the following. 

Considering the five initial values /s, /i^3, /i, fi^2, /i,2,4 on the vertices of the 4-cube (see Fig. ISJ 
one can find the values on all other vertices in a unique way using the dpKdV equation on each 
2-dimensional face because of its 4-dimensional consistency. From these five initial values we form 
the differences 

X = fl- f3, y = h,2- h,3 Z = /i,2,4 - /l,2,3 , (48) 

which are assigned on (the diagonals of) the 2-dimensional faces of the 4-cube (Fig. We note 
that the value /i,2,3 can be expressed already in terms of /s, /i_3, /i, /i^2 through the 3-dimensional 
consistency of dpKdV on the "inner" cube. Next we apply successively the map R : {x, y) i— > (u, v) 
given by (jUj) on the (a) "inner", "front", "left" and (6) "right", "back", "outer" 3-dimensional 
cubes to obtain the following composite maps 

(a) : (x, y, z) — > (x2, z) — > (3:2,4, y-i, ^-i) — > (2:2,4, y-i,4, Z-1-2) , (49) 

(b) : {x, y, z) — > {x, y^, ^-2) — > (2^4, ^4, ^-1,-2) — > (2^2,4, y-i,4, Z-1-2) ■ (50) 

The fact that the two ways of obtaining the values for (2:2,4, y-i, 4, 2-1,-2) lead to identical results, 
and thus to the YB property of the map (jUj), is guaranteed by the 4D consistency. One can notice 
that the evolution of the Yang Baxter variables takes place on two parallel layers of the lattice 
i.e. {(ni,n2,0) , (ni,n2, 1) ,ni,n2 S Z}. This is reflected to the fact that there are six out of the 
eight 3-dimensional faces (cubes) of the 4-cube involved in the compatibility. 

A final comment about the map H47|) is that under the transformation x 1— > —2;"^, v 1— > —v~^, 
u 1-^ Piu, y I— > P2y it becomes the first YB map in Tabled for 6 = 0. Thus the Fjjj map in the 
classification of is also retrieved from dpKdV in Z^, by using the invariants of the symmetry 
subgroup H. 
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3.2 Consistency of dpKdV around a 4-cube and the Harrison map 

We are going to show that the Harrison map appears in a similar manner as previously using 
now the invariants of the full symmetry group G of dpKdV and extending both the equation and 
its symmetry group in Z^. 

The symmetry group G can be naturally extended to the corresponding system in Z^. Now let 
us consider the following invariants of this group: 

fl-fs /l - /4 ( fl- fi\ (fl-fz\ ^g-j^-j 



/2 — /s ' /2 — /4 ' V /2 — /4 / 3 ' \ /2 — /: 

where the subscript i means the shift in the i-th direction. The natural place for them to live in 
are the corresponding 2-dimensional faces. Next we derive the relations between these variables. 

First of all equations (|15j) and their forward shifts with respect to the lattice directions 3 and 
4 imply that the following relations 

(/l - /3)(/l - /4)3 = (/l - /4)(/l - /3)4 , (52a) 

(/2 - /3)(/2 - /4)3 = (/2 - /4)(/2 " /3)4 , (52b) 

(/l,3 - /)(/l,4 - /)3 = (/l,4 - /)(/l,3 - /)4 , (52c) 

(/2,3 - /)(/2,4 - /)3 = (/2,4 " /)(/2,3 " /)4 , (52d) 

hold on the "back" and "left" 3-cubes of the 4-cube depicted in Fig. |21 Dividing memberwise 
equations (|52aj) . (|52b|) and rearranging terms we get 



fi — h [ fi- 











h h 




1/2 


-h) 



(53) 



f2 — fa \f2 — fi. 
In terms of the variables (|5H) equation (|53|) reads 

xy = uv . (54) 

On the other hand, we can rewrite equation (|52c|) in the equivalent form 



which by using equations (fT5|) reads 

ai - "3 Q2 - «3 / /l - /3 \ \ ^ «1 - «4 A _ a2 - Q4 / /l - fi \ j ^gg^J 



/i - /3 V ai-a3 \f2-f3J4J fi- fi V oii-a^ \f2-f 

Multiplying both terms of equation with (ai — 02) (/i — /2) and rearranging terms, the latter 
takes the form 

^ _ 02 - «4 \ A _ /2 - /3 \ A 02-03 ( f\ - fs 



ai-a^J V fi-fsj V ai - ^3 V/2-/3/4 

^^^^ 

ai-asj \ /1-/4/ V ai-a4 V/2-J4/3/ 
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Finally, recalling the defining relations of the variables x, y, u, equation ()57() becomes 

(1 - 72)(1 - - 71 n) = (1 - 7i)(l - - 72 y) , (58) 

where 

71 = , 72 = . (59) 

ai — 03 ai — 04 

A similar calculation starting with ()52d|) and using (|54|) . delivers the same relation (|58j) . 

Proposition 3.1. The invariants h51]) of the symmetry group G of the dpKdV equation extended 
to are related by the YB map, which coincides with the Harrison map 

Indeed, comparing (|54j) and (|5<S|) with the relations from which we obtain the last map in Table 
we deduce that x, y and u, v are related by the Harrison map. The fact that this map satisfies the 
YB relation can be also derived from the consistency property and geometry of the 5-dimensional 
lattice. 

4 Multi-component YB maps 

In this section we show that YB maps can be constructed equally well from 3D consistent multi-field 
discrete equations, for which no classification scheme exploiting the multi-dimensional consistency 
property has been obtained yet. 

4.1 YB map from the discrete modified Boussinesq system 

The discrete modified Boussinesq (dmBSQ) equations [^Tj involve two fields f,g : iP —>■ CP^ and 
are given by the system 

ai /2 - "2 fi g ai fi g2 - 02 /2 gi 

71,2 = g , gi,2 = 7 ■ (60) 

ai gi - a2 52 / «! gi - "2 52 

Its 3D consistency is provided by a lengthy but straightforward calculation that delivers the sym- 
metric values 

<T aiajfkiaigi - ajgj) cr aiajgkiaifj - ajfi) 

Jl,2,3 — / 7 7 7 V' 91,2,3 — g 7 7 7 (,t)J-j 

with respect to any permutation of the indices (1, 2, 3). Here the cyclic sum cr is over the subscripts 

ijk 

{i,j, k) = (1, 2, 3), (3, 1, 2), (2, 3, 1), and similarly the cyclic sum a is over (i, j) = (1, 2), (2, 3), (3, 1). 

ij 

The explicit dependence of /i,2,3j 9i,2,3 on the values f,g implies that dmBSQ does not satisfy the 
so-called tetrahedron property, which is an additional assumpion in the classification scheme in 
for one-field discrete equations. 
Using the symmetry generators 

vi = /5/, V2=gdg, (62) 
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of the dmBSQ equations, we define as YB variables the fohowing joint lattice invariants 

1 /i 1 /l,2 1 /l,2 1 /2 .„„s 

3; = -7' y = u =-^, V = —, (63) 

J Jl J2 J 

2 91 2 51,2 2 51,2 2 52 

X = — , y =—, u = —, V = — . (64) 
5 51 52 5 

It is immediately seen that the above equations imply that 

1111 / \ 

X y = u V , X y = u V . (65) 

Moreover, the lattice equations H6U|) can be expressed in terms of the above invariants as follows 

1 1 aiv"^ — Q2X^ 2 2 aix^v"^ — a2V^x'^ 

U V = ;3 5", u V = 5 :3 — . (66) 

aix^ — a2^^ aix^ — 02^ 



Finally, solving equations H65|) . ()66|) for (m*,w*) we obtain the reversible YB map 

91 9191 99 

L X + a2 x X y + ai Q2 X y 
ai 02 + ai^ x^ x^ yi + 02^ x^ y2 



1 i< 1 1 4-1 A ai a:" + 02" X" X" + ai Q2 X" 
u=yA, u=Xj4 j4 = 5— — „ 2 9 ' ("'^) 



2 2d 220-1 D ai^ a;^ + a2^ a;^ + ai a2 x^ y2 

u =y B, V =x B B = 1 2 1 > 2 2 2 ■ (^'^^ 

02 x^ + ai a2 x^ x^ y^ + ai^ x^ 

4.2 YB map from the discrete potential Boussinesq system 

Discrete potential Boussinesq (dpBSQ) equations is the second member in the so-called lattice 
Gel'fand-Dikii hierarchy PHI- The dpBSQ equations, in the form they were studied recently in 
involve three fields f,g,h : 1? — > CP^, and they are given by the following system 

hi = ffi-g, (68a) 
h2 = ff2-9, (68b) 
h = fh,2-gi,2-^^^. (68c) 

Jl - J2 

For the purposes of the present discussion, equations (|HH|l exhibit the interesting feature that the 
joint invariants of two symmetry generators are enough to construct a YB map. In connection with 
this issue we note that for an elementary Cauchy problem on a staircase, we should impose initial 
values {f,g,h), {fi,gi), (/2,92)> only. From these data the values (/ii,/i2) and (/i,2i 5i,2i ^1,2) are 
determined uniquely. In particular, equations ()68a() . i)68b|) imply that 



f 51-52 .„„N 

/i,2 = 7 7- , (69) 

Jl - h 



and subsequently the values hi^2 and gi^2 are determined from equation H68a|) (or equivalently ()68b|) ) 
and (|68cj) . respectively. 

Using the infinitesimal invariance criterion H33|) to determine the symmetries of equations ()68(l . 
we find that two particular symmetry generators are given by the following vector fields 

^i = df + fdg + fdh, ^2 = dg-dh. (70) 
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They generate the symmetry transformations 

2 2 

■.if,g,h)^{f + eug + eif + ^,h + sif + ^) , (71) 
:{f,g,h)^{f,g + e2,h-e2), (72) 

respectively. We define now as YB variables the following invariants 

= fi- f , = fi,2 - fi , 

= 91-9- /(/i - /) , = 51,2 -91- /i(/i,2 - /: 



1) > 

x^ = hi-h- /(/i - /) , y'^ = /ii,2 - 51 - /i(/i,2 - fi) 



= fl,2 - /2 , = f2- f , 

= 91,2 -92- /2(/l,2 - /2) , v"^ = 92- 9 - /(/2 " /) , 
«3 = _h^_ f^^f^^^ _f^)^ y^ = h2-h- /(/2 - /) . 



(73) 



They are functionally related by 



+ = + y\ (74a) 
+ = + + x^y^ - u^v^ , (74b) 
j^y^ = + x^y^ - u^v^ . (74c) 



X =-y + r T r , (75a) 



Moreover, the system of equations formed by H68ap - ()68c|) . (|68b|) - ()68c() and 1)69^ can be written in 
terms of the above invariants as follows 

x^ — 

x3^^2^y2^^1yl^^lZL^^ (75^^) 
^3^^2^^2^^V^^1ZL^^ (75^) 

x^ — 

Solving the system of equations (|74|) . (|75|) for (n*,f*) we obtain the following YB map 

v)- = y^ — (ai — Q2)r~^ , f ^ = X"^ + («! — Q2)r^-'^ , 

u'^ = y"^ + {ai - a2)[ai - a2 - 2y^T)T~^ , v"^ = x"^ + {ai - a2){x^ + y^)T-^ , (76) 

= + (ai — 02) (ai — 02 + (a:"*^ - y^)T)T''^ , = x^ , 



where F = — + x^y^ + y^. 

4.3 YB map from discrete Calapso equation and nonlinear a— model. 

In a recent study on discrete isothermic surfaces, Schief j2U| introduced the following vector gener- 
alization of the dpKdV 

(/l,2-/)= |"'""|2 (/l-/2), (77) 
1/1 ~ /2I 
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f : ^ C", under the name discrete Calapso equation. Equation (|77() is three dimensional 
consistent since for an initial value configuration /, /j as in Fig. |2fa), one finds that the value 
f 1,2,3 is given by 

where 

A = (qi — a2){ai - Q3) , /i = (ai - a2)(a2 — 03) , = (ai - a3)(a2 — 03) . (79) 

The consistency property is readily checked since /i 2 3 is invariant under any permutation of the 
indices (1,2,3) labeling the field variables and the corresponding lattice parameters. 

The aim now is to construct a YB map from the equation ()77() . Using the translational invariance 
of equation (|77|) . we define the following YB variables 

X = fi-f, y = /l,2-/l, ^t = /l,2-/2> V = f2-f, (80) 

on the edges of a square, which are related by 

X + y = u + V. (81) 
On the other hand, equation (|lflj) can be written in terms of the variables (|8fl|l in the form 

\x -\- y\ 

Hence, equations ()81() . H82|) deliver the following reversible YB map 

^t = y + l — ^ — i^{x + y), v = x-- — ■ — -K{x + y). (83) 
\x + y\^ + 

Moreover, in PD] it was shown that discrete Calapso equation (|77j) can be specialized to an integrable 
discrete version of the 0(n + 2) nonlinear a- model. This reduction is accomplished by imposing 
the constraint 

l/l' = l, (84) 

on the discrete Calapso equation ((77)) . Since the shifted values of / with respect to any lattice 
directions should also satisfy constraint (jHU, equation (|77|) is compatible with this constraint 
whenever 

2/ • /2 - 2/ • /i = ai - a2 . (85) 
This requirement can be satisfied by taking 

-2/ • /i = ai , -2/ • /2 = a2 . (86) 

In terms of the variables (|8U|) . the above constraints translate to 

\x\^ = 2 + ai, \y\^ = 2 + a2. (87) 
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In view of the previous relations the map ()83() obtains the form 



u = y + 



- \y 

\x + 



V = X — 



- \y 

\x + 



{x + y). 



(88) 



By straightforward calculations one finds that 



u 



2 



= \X 



2 



|2 I |2 

v\- = \yy . 



(89) 



Using the above identity it is easily established that the map R : {x, y) (u, v) given by (|88|) 
is a reversible YB map. Up to a permutation this map was first considered by Adler j^H] in the 
geometric problem about recuttings of the polygons. 

5 Perspectives 

We have shown how the symmetry analysis of integrable equations on quad-graphs can be used 
in order to construct YB maps. In particular, we derived the Harrison map from the consistently 
extended discrete potential Korteweg - de Vries equation to the 4-dimensional lattice. The main 
question now is how far this example can be generalized. In particular, for a given multi-parametric 
symmetry group, is there a general relation between the structure of the invariants and the geometry 
of the objects which the YB variables are assigned to ? What is the dimension of the lattice in 
which the discrete equation should be extended to ? The analysis of other equations from [14] may 
clarify these issues. 
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